Gaussian model of explosive percolation in three and higher dimensions 
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The Gaussian model of discontinuous percolation, recently introduced by Araujo and Herrmann 
[Phys. Rev. Lett., 105, 035701 (2010)], is numerically investigated in three dimensions, disclosing a 
discontinuous transition. For the simple-cubic lattice, in the thermodynamic limit, we report a finite 
jump of the order parameter, J = 0.415 ± 0.005. The largest cluster at the threshold is compact, 
but its external perimeter is fractal with fractal dimension dA = 2.5 di 0.2. The study is extended 
to hypercubic lattices up to six dimensions and to the mean- field limit (infinite dimension). We 
find that, in all considered dimensions, the percolation transition is discontinuous. The value of the 
jump in the order parameter, the maximum of the second moment, and the percolation threshold are 
analyzed, revealing interesting features of the transition and corroborating its discontinuous nature 
in all considered dimensions. We also show that the fractal dimension of the external perimeter, for 
any dimension, is consistent with the one from bridge percolation and establish a lower bound for 
the percolation threshold of discontinuous models with finite number of clusters at the threshold. 

PACS numbers: 64. 60. ah, 64.60.al, 89. 75. Da 



I. INTRODUCTION 

Percolation, one of the most famous models in sta- 
tistical physics, has been extensively considered as a 
paradigm to study connectivity and transport [l]. Re- 
cently, Achlioptas, D'Souza, and Spencer [2] have pro- 
posed a best-of-two product rule for bond selection char- 
acterized by a more pronounced transition than in the 
random case, being apparently discontinuous [3]. This 
model has been analyzed on several different graphs [11 15] 
and the ambiguous reported results raised controversy 
about the nature of the transition j6HT2] , with analytical 
[6l[T2^ and numerical [5 , 6 results showing the continuous 
nature of the transition in the original best-of-two prod- 
uct rule. Several different models have been studied to 
shed light on the main mechanisms leading to a discon- 
tinuous percolation transition p!3UT6) . A generalization 
to a best-of-m product rule has also been proposed [TT] 
and a tricritical point found when explosive percolation, 
obtained with m = 10, is diluted with classical percola- 
tion on a square lattice [18]. 

Araujo and Herrmann [16] introduced two models 
yielding clear discontinuous transitions: the largest clus- 
ter and the Gaussian models. The study of the former 
discloses the control of the largest cluster as a way to 
obtain homogenization of the cluster sizes and, conse- 
quently, an abrupt transition. Since the properties of the 
best-of-two product rule depend crucially on the topol- 
ogy [Ull] 5 in this work, we study the Gaussian model on 
hypercubic lattices up to dimension six and in the mean- 
field limit (infinite dimension). We report that, for all 
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dimensions, the Gaussian rule leads to a discontinuous 
transition at the percolation threshold and that the frac- 
tal dimension of the largest-cluster external perimeter is 
compatible with the one reported for bridge percolation 

This manuscript is organized in the following way. In 
the next section we describe the Gaussian model and an- 
alyze its properties on the simple-cubic lattice. In Sec. [Hi 
the study of the model is extended to higher dimensions. 
We leave the final remarks for Sec. lIVI 



II. THE GAUSSIAN MODEL ON THE 
SIMPLE-CUBIC LATTICE 

We start by considering a simple-cubic lattice with lin- 
ear size L and periodic boundary conditions in all direc- 
tions. In the initial configuration, all the 3A^ bonds are 
empty, such that there are N = clusters of size unity. 
At each iteration, a new bond is randomly chosen among 
the empty bonds and occupied with probability 
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where s is the size of the cluster that would be formed by 
occupying the selected bond and s is the average number 
of sites per cluster if the bond would be occupied. For 
bonds which connect sites belonging to the same cluster, 
s is taken as twice the size of the cluster, a is a param- 
eter of the model which, for the sake of simplicity, we 
take equal to unity. The proposed method promotes the 
homogenization of the cluster sizes by suppressing the 
formation of clusters differing significantly, in size, from 
the average. 

The difference between classical percolation and the 
Gaussian model can be seen qualitatively in Fig.[T] where 
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FIG. 1. (color online) Snapshots of (a) classical percolation 
and (b) the Gaussian model of discontinuous percolation, at 
the percolation threshold, on a simple-cubic lattice with 128*^ 
sites. To enhance visibility, the front cube of size 64^ has 
been left out in both pictures. While in the classical case 
clusters have very different size and a fractal shape, for the 
Gaussian model clusters are rather compact and with a char- 
acteristic size. The seven largest clusters of the configurations 
are shown in (c) for classical percolation and in (d) for the 
Gaussian model. 



we show snapshots for both models, on the simple-cubic 
lattice, at the respective percolation thresholds. For clas- 
sical percolation, Fig.jlja) and (c), the clusters are fractal 
and of very different sizes, following a power-law distri- 
bution [l], whereas for the Gaussian model (Fig.jljb) and 
(d)) clusters are rather compact and of comparable size. 
Note that, while for classical percolation, the number 
of clusters is large, ^ 0.27A^, the number of clusters in 
the Gaussian model is significantly smaller. In fact, as 
we show here, in the thermodynamic limit, the Gaussian 
model is characterized by a finite number of macroscopic 
clusters at the threshold. 

To analyze the behavior of the order parameter, 
namely the fraction of sites in the largest cluster, we 
measure, for each sample, its jump J, defined as the 
maximum change obtained as one sequentially occupies 
bonds in the system [9^, T3^ . For every considered linear 
system size L, we average the jump and the fraction of 
occupied bonds p, at which it occurs Pc,j over several 
configurations. We take the latter as an estimator for 
the threshold in the thermodynamic limit. Recently, Lee 
et al. [To] defined it as the upper pseudo-transition point 
and used it to pin down the threshold. Plotting J as a 
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FIG. 2. (color online) System size dependence of the maxi- 
mum jump of the order parameter, J (■), for the Gaussian 
model of discontinuous percolation on the simple-cubic lat- 
tice. In the classical case, the change in the order parameter 
shrinks with the system size and is zero in the thermodynamic 
limit (not shown), whereas in the discontinuous case, a finite, 
non-zero value (J = 0.415 =b 0.005) is obtained in this limit. 
Note that, whereas for the product rule and classical perco- 
lation the size of the jump decreases with the system size, 
for the Gaussian model it even slightly increases for the same 
range of system sizes. Results have been averaged over 10^ 
samples for the smallest system size and 1.1 x 10^ samples for 
the largest one. Random numbers have been generated with 
the algorithm proposed in Ref. [20]. To identify the clusters 
and keep track of their properties we have considered the la- 
beling scheme proposed by Newman and Ziff [21], related to 
the Hoshen-Kopelman algorithm [22] . 



function of L reveals that for the Gaussian model, in 
the thermodynamic limit, the jump has a finite value of 
J = 0.415 ± 0.005 (see Fig.[2|, as expected for a discon- 
tinuous transition. This result is in contrast to the ones 
for classical percolation and the product rule where, for 
the same range of system sizes, the size of the jump di- 
minishes and eventually vanishes in the thermodynamic 
limit [9]. 

To determine the threshold two different estimators 
have been considered: the average fraction of occupied 
bonds at which the jump occurs Pc,J [S [E], and the 
position pc,M [23^, of the maximum in the second moment 
of the cluster size distribution, excluding the contribution 
of the largest cluster (of size Smax), 

M^ = M2-sLx/iV , (2) 

where M2 = ^iS^/N and Si is the size of cluster i. 
Figure |3] shows the system size dependence of both es- 
timators on the simple-cubic lattice. Asymptotically, a 
dependence on is found, compatible with using the 
same exponent a — 1.69 ± 0.10 for pc,j and Pc,M- The 
estimators are extrapolated to the thermodynamic limit, 
given by 0, and combining both methods yields 

Pc = 0.3468 ± 0.0005. Note that, for an equilibrium first- 
order transition, a = d [24|. To shed light on the obtained 
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FIG. 3. (color online) Threshold pc, as a function of the 
inverse linear system size for the Gaussian model of dis- 
continuous percolation on the simple-cubic lattice. pc,j (x) 
stands for the average fraction of occupied bonds at which the 
jump occurs and Pc,M (+) for the position of the maximum 
of the second moment of the cluster size distribution with- 
out the contribution of the largest cluster. The percolation 
threshold is estimated to be 0.3468 =b 0.0005. Results have 
been averaged over 10^ samples for the smallest system size 
(32^ sites) and 10^ samples for the largest one (256^ sites). 
The inset shows the same for the case where only merging 
bonds are considered; this is, all clusters are trees (loopless). 
In this case, the threshold estimators [pc,j (•), Pc,m (a)] de- 
pend asymptotically linearly on . The threshold for the 
loopless case is estimated to be 0.3333 ± 0.0004. Results have 
been averages over 10^ samples for the smallest system site 
(50^ sites) and 10^ samples for the largest one (512^ sites). 
Error bars are smaller than the symbol size. 



value, we investigate the dependence of the threshold es- 
timators on the system size, under the constraint that 
only merging bonds are considered; this is, all clusters 
are trees (loopless). In the inset of Fig. [s] we see Pc,j and 
Pc,M as a function of . One observes that in this case, 
a = d. Therefore, in the Gaussian model a differs from d 
due to internal bonds which do not influence the cluster 
structure or the size of the jump. 

In Fig.[4] we see the size dependence of the maximum of 
the second moment per lattice site. For every sample, we 
measure the maximum M2{Pc,m) / and average over all 
samples. For large system sizes, this quantity is constant, 
as expected for a discontinuous transition. 

The scaling behavior of the standard deviation of the 
order parameter, defined as 



FIG. 4. (color online) Maximum of the second moment of 
the cluster size distribution per lattice site, M2{Pc,m) / (*), 
as a function of the linear system size L. The second moment 
per lattice site tends toward a constant value, as expected for 
a discontinuous transition. Results have been averaged over 
10^ samples for the smallest system size and 2.4 x 10^ samples 
for the largest one. 
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FIG. 5. (color online) Standard deviation of the order pa- 
rameter Xoo as a function of the bond occupation fraction p 
for different linear system sizes L. One observes that the peak 
increases and narrows with the system size. In the inset we 
see Xoo as a function of the scaling variable {p — pc)L"' ^ with 
a = 1.64, for different linear system sizes L\ 16 (+), 32 (x), 
64 (*), 128 (□), 256 (•), and 512 (■). Results have been av- 
eraged over 10^ samples for the smallest system size and 28 
samples for the largest one. 



is shown in Fig. [5] For large systems, the maximum of 
Xoo tends toward a constant value. These results are 
a strong evidence of a discontinuous transition since a 
nonzero value of Xoo is obtained at the transition point, 
as expected in the presence of a jump in the order pa- 
rameter [24l [25] . In addition, the plot is consistent with 
the exponent a = 1.69 ±0.10 and = 0.3468 ±0.0005, as 



determined from the finite size scaling of Pc, J and Pc,M- 

We measure at Pc the external perimeter of the largest 
cluster A. The external perimeter is defined as the num- 
ber of sites which do not belong to the largest cluster but 
are nearest neighbors of sites in this cluster [26] . One ob- 
serves that, at the threshold, the external perimeter of 
the largest cluster scales asymptotically with the system 
size as A ~ L^^, where dA = 2.5 ± 0.2 (see also Fig.js]). 
On the square lattice, the fractal dimension of the exter- 
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nal perimeter was shown to be related to several other 
models [161 [H] • The value reported here for the simple- 
cubic lattice agrees within its error bars with the one for 
watersheds and the optimal path cracking |28j as well as 
with the set of bridges in bridge percolation [19] . Clus- 
ters at the threshold are compact with fractal external 
perimeter, as was also reported for 2D [16 and for irre- 
versible aggregation at high concentration [29] . 



III. HIGHER DIMENSIONS AND MEAN-FIELD 
BEHAVIOR 

The Gaussian model yields a discontinuous percola- 
tion transition in two and three dimensions. How does 
the nature of the transition depend on the dimensionality 
of the system? To address this question, we consider the 
Gaussian model on hypercubic lattices up to d = 6, the 
upper critical dimension of classical percolation [30j. In 
addition, the mean-field behavior of the Gaussian model 
is investigated. In the latter case, we take a system with 
N sites which can be fully interconnected giving a total 
of N{N — l)/2 links, and we add links between sites with 
probability given by Eq. ([T]) . For this system, p is de- 
fined as the average number of links per site. Occupying 
links randomly, without any additional rule, would re- 
cover Erdos-Renyi percolation, where Pc = 1/2 (see, for 
example, Ref. [31 ). 

Figure [6] shows the jump J, as a function of the inverse 
system size A/""^, for 3 to 6 dimensions and for mean- field. 
We observe that, in the thermodynamic limit, J has, 
within the error bars, the same finite value in all consid- 
ered dimensions, consistent with the value found in three 
dimensions, J = 0.415 ±0.005 (see Fig.|2|. In general, we 
expect for a discontinuous percolation transition to find 
few macroscopic clusters at the threshold, as initially dis- 
cussed by Friedman and Landsberg |3]. Nagler, Levina, 
and Timme [9 have added that, for strongly discontin- 
uous transitions, where the largest cluster cannot grow 
directly, the number of clusters is finite and the transi- 
tion occurs when the two largest clusters merge. The 
jump is then bounded by two limits: either the clusters 
have the same size, giving J = 1/2, which is the largest 
possible jump size in discontinuous percolation, or the 
largest cluster is of size ^2/3 and the second largest of 
size ~ 1/3, giving J = 1/3. The latter case corresponds 
to situations where the second cluster merges with the 
third one and becomes the largest one, of size 2/3, called 
overtaking in Ref. [9 , merging later with the one of 1/3. 
On the other hand, the former results from four clusters 
of equal size which merge in pairs. This case is expected 
for the global competition proposed in Ref. [13 , in the 
mean- field limit. The Gaussian model at any dimension 
also promotes the homogenization of the cluster sizes and 
the values of the jump are within the proposed interval. 
The same idea can be considered to understand the be- 
havior of the maximum of M2/N ^ taking place at Pc,M 
which is our second estimator. 
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FIG. 6. (color online) System size dependence of the jump 
J, for the Gaussian model of discontinuous percolation on the 
hypercubic lattice of dimension 3 (x), 4 (♦), 5 (a), and 6 (•), 
as well as in the mean-field case (*). In the limit 0, 
the jump has within the error bars the same finite value 0.415 
for all considered graphs. The solid line is a guide to the eye 
and of the form 0.415 — AQN~^ . For the sake of comparison, 
we plot the jump as a function of the inverse system size . 
Results have been averaged over 10^ samples for the smallest 
system size and at least 10 samples for the largest one. 
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FIG. 7. (color online) Cluster-size distribution for the Gaus- 
sian model on the simple-cubic lattice. The fraction Us (•) of 
clusters of size s times the size of the cluster s/N \s shown 
as a function of s/N. One observes that the distribution is 
bimodal, as expected for discontinuous transitions. The sys- 
tem size is 64^ sites, results have been averaged over 1.6 x 10^ 
samples, and error bars are indicated. 



In Fig.[7| we see the cluster-size distribution for the 
Gaussian model on the simple-cubic lattice at the perco- 
lation threshold, p = p^ As previously observed in 2D 
[16] . a bimodal distribution is obtained, in contrast with 
the power-law behavior observed for random percolation 
[1 and the best-of-two product rule [4]. Since the contri- 
bution of the largest cluster is neglected, there is a cut-off 
at s/N = 0.5. 
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As in three dimensions, we also determine the perco- 
lation threshold for the Gaussian model in dimensions 4, 
5, and 6, as well as in the mean- field limit, by combin- 
ing both estimators: Pc,j and Pc,M- Table [l| shows the 
threshold values for different dimensions. One ob- 
serves that Pc decreases with the dimension, though it 
remains always above the values for classical percolation 
[32 (shown in the same table for comparison). For the 
Gaussian model in the mean-field limit we find Pc to be 
compatible with unity but note that, in this case, p is 
defined as the fraction of bonds per site and not the frac- 
tion of occupied bonds, as in the lattice case. Below we 
establish a lower bound for the Pc of models yielding a 
discontinuous percolation transition with finite number 
of clusters at the threshold. 

Consider an arbitrary percolation model which starts 
with isolated clusters of unit size, adding bonds se- 
quentially until a certain fraction of occupied bonds is 
reached. Let us denote by c{p) the number of clusters at 
a given fraction of occupied bonds p. At each iteration, 
added bonds to the system can be merging bonds - con- 
necting two clusters - or redundant bonds - connecting 
nodes of the same cluster Only the former bonds 

change c{p). The number of clusters reduces by one if 
the bond is a merging bond and does not change if it is 
a redundant bond. If we define r{p) as the probability 
that an added bond is redundant, then 



db ^ 



(4) 



where b = pNd is the number of occupied bonds in a 
(i-dimensional hypercubic lattice. We now take the limit 



lim 

N^oc 



dc 
dpN 



= -d[l-r{p)] 



Integrating over the interval < p < Pc gives 



N 



lim [c{pc) - c(0)] /N = -d / dp [1 - r{p)] 



(5) 



(6) 



TABLE L Percolation threshold pc for the Gaussian model 
of discontinuous percolation on the hypercubic lattice of di- 
mension d and in the mean- field limit. For comparison, the 
percolation thresholds for classical percolation are shown in 
the third column. Note that, in this table, for all models, p 
is defined as the fraction of occupied bonds, p = t/{Nd) for 
hypercubic lattices and p — t/N for mean field, where t is the 
number of occupied bonds in the system. 
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FIG. 8. (color online) System size dependence of the exter- 
nal perimeter of the largest cluster A, at the threshold, for 
the Gaussian model of discontinuous percolation on the hy- 
percubic lattice of dimension 3 (*), 4 (♦), 5 (•), and 6 (■). 
One observes that A asymptotically scales with the system 
size as A ^ L^^. The solid lines have slopes of 5.9 it 0.8, 
4.9 ± 0.7, 3.6 =b 0.4, and 2.5 ± 0.2, respectively Results have 
been averaged over 10^ samples. 




This equation is valid for any percolation model regard- 
less the nature of the transition. For example, the per- 
colation threshold for the classical tree-like case can be 
obtained by taking r(p) = and the proper number of 
clusters at the threshold [33ll35] . 

Assuming the cases where the homogenization of the 
cluster sizes leads to a finite number of clusters at jl] , 
and since c(0) = A/", 



Pc l/d- 



r 

Jo 



dpr{p) 



(7) 



we obtain that Pc > l/d. Note that for tree-like models 
r{p) and, if c{pc)/N as N ^ oo, Pc = l/d. 
For the Hamiltonian model of explosive percolation, in- 
troduced by Moreira et al. [14 , the same result was 
derived in the mean-field limit in an independent way 



and numerically observed in the lattice. The value re- 
ported by Manna and Chatterjee [13] for the case with 
global competition is also consistent with this result. 
Both the lower bound for the threshold and the solu- 
tion c{p < Pc)/N ^ 1 — pd, obtained for vanishing small 
probability of redundant bonds, are consistent with the 
numerical results for the Gaussian model. For increasing 
dimension the contribution of redundant bonds decreases 
and Pc approaches l/d (compare Tab.[l|. In the mean- 
field limit, this asymptotic behavior also agrees within 
error bars with the reported results. The transition is ob- 
tained when the number of added bonds equals the num- 
ber of nodes N . Since the maximum number of bonds is 
N{N — l)/2, the fraction is zero in the thermodynamic 
limit. 

Measuring the size dependence of the largest-cluster 
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external perimeter A, at the percolation threshold its 
fractal dimension is obtained for dimensions 4, 5, and 
6 (see Fig.|8|. For increasing dimension, dA seems to 
approach d (see Tab.|l|. These exponents are compatible 
with the ones found for bridge percolation, corroborating 
their equivalence [19 . 

IV. FINAL REMARKS 

In summary, in this work we studied the Gaussian 
model of discontinuous percolation in three and higher 
dimensions. We disclose that, for any considered dimen- 
sion, the percolation transition is abrupt and character- 
ized by a discontinuity in the order parameter, which 
within error bars is independent on dimension. We iden- 
tify the homogenization of cluster sizes and favoring of 
merging bonds as the key mechanisms leading to such 
an abrupt transition [9l [131 EH Ifi.- For discontinuous 
percolation models with a finite number of macroscopic 
clusters at the threshold, we establish a lower bound for 
Pc as well as a relation between pc and the probability of 
selecting a redundant bond. Studying different dimen- 
sions we show that clusters are compact with a fractal 
perimeter with the same dimension as bridge percolation 
[19], which is also related to watersheds and the optimal 
path cracking [28 . Although all numerical indications 
point in that direction we have no formal proof whether 
the upper-critical dimension for the Gaussian model is 



six, like in the classical case. In addition, the meaning of 
the non-trivial finite-size scaling exponent a = 1.69±0.10, 
consistent for both estimators of Pc^ is still puzzling; an 
analytical treatment of this exponent would be interest- 
ing. Studies of this model have taken a = 1, in Eq. ([T]), 
but it would be interesting to investigate other cases since 
for a = the model boils down to the classical percola- 
tion model. It would be interesting to study how the de- 
scribed properties depend on a. Future work might also 
consist in studying the behavior of other models of ex- 
plosive percolation, like the global competition proposed 
by Manna and Chatterjee [13] and the BFW model dis- 
cussed by Chen and D'Souza pT, in different dimensions. 
Besides, all known models with a discontinuous transition 
imply global information. It is still an open question if 
a discontinuous percolation transition can be obtained 
with only local rules. 
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